Understanding ferromagnetism in Cr-based 3d-5d double perovskites 
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Ferromagnetism in Cr-based double perovskites is analyzed using effective model as well as sim- 
ulation approach. Starting from a microscopic model proposed recently for this class of double 
perovskites, an effective spin-only model is derived in the limit of large exchange coupling at the 
B-site. Analytic expressions for the resultant exchange coupling is derived in certain limiting cases. 
The behaviour of this exchange is used to provide a plausible explanation for the enhanced ferro- 
magnetic tendency as also the enigmatic increase in Curie temperature observed in Cr-based DP-s, 
in the series Sr2CrW06,Sr2CrRe06 to Sr2CrOs06. The superexchange between neighbouring B and 
B' sites is found to play a crucial role both in stabilizing ferromagnetism, especially in the latter 
two compounds, as well as increasing the Tc. 

PACS numbers: 



I. INTRODUCTION 



In recent times, ferromagnetic transition metal com- 
pounds with high Tc have come into the focus for 
their importance in spintronics and other technologi- 
cal applications. While doped rare earth manganites 
have been studied for decades for their CMR property, 
no less well known are the double perovskites, with a 
general formula A2BB'06, A=alkaline/rare-earth metals 
and B/B'=transition metals. The most well known mem- 
ber, Sr2FeMo06 has a Tc of 410K, higher than most man- 



ganitesii^. Moreover, the half-metallic nature^ii, cou- 
pled with the substantial tunnelling magnetoresistance 
obtained in these compounds at low temperatures, espe- 
cially when powdered, makes it valuable for spintronic 
applications^"—. The high Tc implies a strong polariza- 
tion even at room temperatures, enhancing its industrial 
importance. In an effort to boost the Tc even further, re- 
searchers have electron doped this compound^"— li^. Al- 
though the Tc does increase in this process upto some 
point^"'^^, but it has been shown in a recent work that 
upon overdoping, the Tc actually decreases, and the fer- 
romagnetism becomes unstablei^ii^! In fact, it gets re- 
placed by antiferromagnetic phases. This conclusion had 
been derived both within model Hamiltonian approach^^ 
as well as within abinitio approach. Within the model 
approach, it was observed that the effective exchange in- 
teraction between Fe S — 5/2 core spins changes sign 
as the filling increases, signalling a crossover from ferro 
to antiferro. Within the abinitio approach, considering 
the Sr2_2:La2:FeMo06 series^^, it was observed that the 
ferromagnetic state became progressively unstable as the 
number of valence electrons was increased through in- 
creased doping of La, the balance tilting at an electron 
count of about 2.4. While the La-overdoped regime of 
this compound has not been investigated in detail as yet, 
recent experimental data support our claims regarding 
the disappearance of ferromagnetism on electron dop- 
ingi^. Hence simple electron doping by A-site cataion 
substitution is not a very promising method for increas- 
ing Tc beyond a point. There is however, another well 



known technique to change the electron filling, namely to 
substitute instead the B' site ion. In this case, it is well 
known that in Cr-based double perovskites: Sr2CrB'06, 
as one goes across the period, substituting B' in turn 
as W, Re and Cr, the Tc increases progressively, reach- 
ing a high of over TOOK for the Os compound. Obvi- 
ously, the ferromagnetism is not becoming unstable in 
this series of compounds even though the filling is in- 
creasing from 1 to 2 to 3 as one goes from W to Re to 
Os respectively. On the contrary, it is becoming more 
stabilized, as signified by the increasing Tc- In this pa- 
per, we probe the role of a novel superexchange mecha- 
nism, in addition to the kinetic energy driven mechanism 
already prevalent, to account for this anomalous stabi- 
lization of ferromagnetic behaviour. While abinitio and 
variational approaches had earlier-i^ provided a pointer 
to this superexchange mechanism, in this communica- 
tion, we conclusively demonstrate its importance using 
effective exchange calculations and direct numerical sim- 
ulations. The organization of the paper is as follows. In 
the next section, we shall briefly summarize the Hamil- 
tonian as well as main results for the Sr2-xLa2:FeMo06 
series (i.e., cataion site substitution), just to serve as a 
reminder, and for constrast with the Cr-based series (B- 
site substitution). Next, we provide a brief account of the 
main results of abinitio studies upon these compounds, 
which have been reported in detail earlier. In the fourth 
section, we motivate a new modified Hamiltonian for this 
Cr-series. We proceed to derive a low-energy, spin-only 
model from this fermionic Hamiltonian in the next sec- 
tion, and analyze its magnetic properties. 



II. A-SITE CATAION DOPING: 
SUMMARY 



BRIEF 



A-site cataion substitution can be done upon the par- 
ent compound Sr2FeMo06 using, for example La which 
has a nominal valence state of 3-1- , in place of Sr which 
has a valence state of 2+. This would correspond to 
electron doping of the system. Abinitio studies on this 
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series of compoundalA, using NMTO^*' downfolding have 
shown that the relative positions of Fe and Mo or- 
bitals remain almost unchanged upon La doping, only 
the Fermi energy shifts almost like a rigid band picture. 
Thereupon, total energy calculations using Vienna Abi- 
nitio Simulation Package (VASP)^ showed that the mag- 
netic ground state changed from ferro to antiferro as the 
electron doping increased. The spin splitting at the Fe 
and Mo site was also obtained by downfolding everything 
except the t2g orbitals. It was observed that the spin- 
splitting at the Mo-site increased proportional to the fill- 
ing, so that the Stoner / remained essentially constant. 
For example, in Sr2FeMo06 which has 1 electron per site 
the spin-splitting at the Mo site is about 0.13 eV, while 
in La2FeMo06, with a filling of 3, the splitting is 0.37 
eV. This showed that the moment at the Mo site is 
completely induced from the Cr spin, and not intrinsic. 
It was shown in that paper that the Hamiltonian which 
captures this behaviour is given by: 

H = epe^ fLafi'^a + (Mo ^ m\^^mi^a 
ieB ieB' 

-tpM ^ fl,a^J<T,a - tMM ^ "^l^,a"^J<T,a 

-tpp ^ fl„^Jja,a + Jl^^i- fla^afshli (1) 
<ij>(T,a ieA 

The /'s refer to the Fe sites and the m's to the Mo 
sites. tpM, tMM, tpp represent the nearest neighbor Fe- 
Mo, second nearest neighbor Mo-Mo and Fe-Fe hoppings 
respectively, the largest hopping being given hy tpM- o 
is the spin index and a is the orbital index that spans 
the Xig manifold. The difference between the ionic levels, 
^ = epe — (Mo, defines the charge transfer energy. Since 
among the crystal-field split d levels of Fe and Mo, only 
the relevant t2g orbitals are retained, giving rise to on- 
site and hopping matrices of dimension 3x3. The Si are 
'classical' (large S) core spins at the B site, coupled to 
the itinerant B electrons through a coupling Ji 3> tpM- 
Variants of this two-sublattice Kondo lattice model has 
been considered by several author o"'^°i^^~ '— in the context 
of double perovskites. 

The abinitio calculations showed that the charge trans- 
fer energy A remains almost fixed upon La-doping. Con- 
sequently, the doping only corresponds to increasing the 
Fermi energy keeping the levels fixed, thereby increas- 
ing the electron filling. If one further integrates out, or 
downfolds, the Mo t2g degrees of freedom represented by 
771^, and traces over the /^-s, then one can arrive at an 
effective model consisting of only the Iron core spins S", 
and calculate the resultant exchange interaction between 
theroi^. This exchange, when plotted at constant A as 
a function of filling alone, exhibits a change of sign as 
the filling increases, as shown in Refl'^. Thus, the ground 
state changes from ferromagnetic to antiferromagnetic at 
a filling of about 2.2-2.4. 



III. B' SITE SUBSTITUTION: MAIN ABINITIO 
RESULTS 

In the series of compounds Sr2CrB'06, the B^ ion 
which is in nominal 5-1- valence state corresponds to 
bd} ,?>(P,b(P configuration of W, Re and Os respectively. 
These materials have been reported to have Tc as high as 
about 450K, 620K and 725K. with a progressive increase 
as one moves from Sr2CrW06 with one valence electron 
to Sr2CrOs06 with three valence electrons. Taking the 
number of valence electrons as sole consideration, the 
situation of Sr2CrW06,Sr2CrRe06,Sr2CrOs06 is com- 
parable to Sr2FeMo06, SrLaFeMoOe and La2FeMo06. 
However, unlike the Cr-B^ {B^ =W,Re,Os) series, for the 
Sr2_a;La'2;FeMo06, the ferromagnetic was found to 
decrease with increasing La concentration (i.e., increas- 
ing number of valence electrons) and finally the antifer- 
romagnetic phase taking over the ferromagnetic phase. 
The Cr-B/ (B^ =W,Re,Os) series though bear two funda- 
mental differences compared to Sr2_a;Laj;FeMo06 series. 
Firstly, the B^ ions in Cr-B^ being 5d transition metals 
exhibit significant spin-orbit coupling which make these 
materials suitable for magneto-optic applications with 
large signal as has been discussed in Refl. Secondly, three 
different chemical elements, namely W, Re and Os are 
involved in Cr-B/ series while for Sr2-a;Laa;FeMo06 se- 
ries the increased electron count is achieved without any 
changes in the B-B^ sublattice. First principle calcula- 
tions using generalized gradient approximation (GGA)^-, 
carried out using (VASP) as well as Linear Muffin Tin 
Orbital (LMTO) have been carried out for the total en- 
ergy calculations. Thereafter, a few-band, low-energy 
tight binding Hamiltonian was derived using N-th or- 
der muflin-tin orbital (NMT0)i2 formalism. Details have 
been reported elsewhere^^; here we recapitulate the main 
results. 

It was observed that the moment on the B' site in- 
creased from 0.3/iB for tungsten (W) to O.Sl/is for Rhe- 
nium (Re) to lAA^B for Osmium. The total moment, on 
the other hand, goes from 2/Lt_B in W to l^s in Re to Oub 
in Os. It is observed that although the total moment de- 
creases in steps commensurate with the filling, and the Cr 
moment remains almost fixed, the moment at the B' site 
increases drastically, pointing to the growing localization 
of electrons on this site. 

To probe this issue in more detail, the Sr, Oxygen or- 
bitals as well as the eg orbitals of Cr and B' ions were 
downfolded using NMTO formalism, keeping only the Cr 
and B' t2g orbitals. It was observed that the spin split- 
ting at the B' site increased from 0.06 eV at the W site 
to 0.31 at the Re site to 0.53 at the Os site. This is 
obviously a much faster increase than the filling would 
dictate, indicating that the Stoner / is itself increasing. 
If we multiply the splitting for W, i.e., 0.06 by 2 and 3 
respectively, then the extra amount must correspond to 
a different energy scale in the problem. Let us call this 
energy scale as J2, while the traditional extremely large 
spin splitting at the B site, or Chromium site defines the 
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other exchange energy scale, Ji. 

Another important point which emerges out of abini- 
tio calculations is the progressive increase in the charge 
transfer energy A as one goes across the series from W to 
Rc to Os. It increases from 0.51 eV in W, to 0.9 in Re, 
to 1.35 in Os. This situation is to be contrasted with the 
nearly constant A in the Sr2-2;Laa;FeMo06 series. Thus, 
it may be summarized that both A and J2 are increasing 
in the Cr-series from W to Re to Os compounds, although 
J2 is usually much less than A. 



IV. HAMILTONIAN 

From the above considerations, it is obvious that there 
exists a different exchange energy scale to the problem, 
apart from the spin splitting at the Cr site Ji, which, 
being related to the Hund coupling, for all practical pur- 
poses can be considered infinite. The increasing localized 
character of the moment on the B' site makes superex- 
change an obvious candidate for this new energy scale. 
Hence we have added a superexchange term to the Hamil- 
tonian corresponding to superexchange between the t2g 
spin on the B site, and the electron spin on the B' site. 
Hence, the representative Hamiltonian is given by: 



ieB 



ieB' 



(Sr2_a;Laa;FeMo06) using the procedure of SCR. In the 
action corresponding to this Hamiltonian, all but the last 
term is translationally invariant, hence can be fourier 
transformed and written in terms of the bare Green's 
function. Thereafter, the Molybdenum degrees of free- 
dom can be integrated out, and the Iron degrees of free- 
dom traced over, after taking the limit of Ji — >■ cxd, to 
obtain an effective model containing only the Iron spins. 
One may imagine that the same procedure may be fol- 
lowed in this case. There is a fundamental difficulty 
in this case though, due to the fact that there is spin- 
disorder even on the Molybdenum site, making it impos- 
sible to diagonalize and integrate out the Molybdenum 
degrees of freedom. Hence, the entire procedure of thin- 
ning of degrees of freedom has to be done in real space. 
The action corresponding to the Hamiltonian in Eqn[2]is 
given by: 



<ij> 



(3) 



B'B/ E 

<ij>(T,a. 



tec Y fL,af3'y,a + Jl ^ S,; • //^Tq/j/, 



<.ij'>(7,a. 



j/3 



ieA 



ieB 



Now, since the spin orbit coupling at the B' site is large 
for the 5d elements like Re and Os compared to the 4d 
elements like Mo, the onsite spin anisotropy is also ex- 
pected to be high. Let us choose our axis of quantization 
along this anisotropy axis, whereupon we need only con- 
sider the diagonal components of the last J2 term. Hence, 
the action can be written as: 



(2) -^ = E 



This spin-fermion Hamiltonian is in general difficult to 
solve exactly, but some headway in the direction of un- 
derstanding the low-energy magnetic behaviour can be 
made if on can obtain a spin-only model from this with 
an effective exchange. In principle this should be possible 
by tracing out the fermion degrees of freedom, but that 
is a herculean task considering the multitude of possi- 
ble configurations. Instead, we consider the approximate 
but enlightening procedure of Self-Consistent Renormal- 
ization (SCR) devised by Kumar and Majumdar— . 



V. DERIVATION OF EFFECTIVE SPIN MODEL 



(4) 



where the matrix M is given by: M — {iuJn ~ cmo)! + 
•hS^j^^cTz- Hence, first taking the Ji — ^ 00 limit^, and 
then integrating the B' degrees of freedom out, we get 
the following form for the action: 



In an earlier paper^'^, I had obtained an ef- ^2 

fective spin model from the two-sublattice Kondo — ^ £M -cos—cos—jj^jj^. 

lattice model(Eqn [T|) appropriate for the Fe-series ^^.^^^^ — ^Mo + J2J25 ^ ^ 
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E ■ 

<<jj>> 



'-FM 



-szn— Sin— 7,' 7j-„ 



(5) 

where 7i_„ refer to transformed spinless Ferraion opera- 
tors at the Iron site, while 9i refer to the azimuthal angle 
made by the Fe core spins with the z-axis. Following 
usual practice, we have neglected the Berry's phase de- 
grees of freedom connected with the polar angles. 

Hence, if we consider total energy U = — ^§^^ as 

befor&i3.^ then we get: 



2 

FM 



<<ij» 



calculate the exchange as follows: 



-4 = ^1 iUn 



'-FM 



A -Jo 



Now we calculate U 



to obtain 



t%j^,i{2iujn - A - J2) 



{iuj„ - A - J2)2 
Evaluating the Matsubara sum. 



(7) 



(8) 



(2zc^„- A- JaE^^I) 



„ stn — sm— 
2 2 2 



(6) 



So far, no approximation has been made except that 
J2 <C Ji and a large anistropy at the B' site. How- 
ever, to actually calculate the exchange from this ex- 
pression, one first needs to start with some spin back- 
ground, say random, and calculate the Iron Green's func- 
tion in real space in this background using the original 
Hamiltonian [5] Then, one needs to recalculate the spin 
background from the effective spin Hamiltonian [B] again 
using some technique like Monte Carlo. This procedure, 
to be repeated till self-consistency, defines the technique 
of Self-consistent Renormalization^^. The entire proce- 
dure is however numerically expensive and difficult. We 
shall, in the next section, instead, try to calculate the ex- 
change analytically in certain simple arrangement of the 
background spins: namely fully ordered ferromagnetic 
configuration. This can act as the first step of the SCR 
procedure, which will be taken up in a future work. As 
we shall see, it will give us an important analytic pointer 
at the actual magnetic ground states possible in these 
compounds. 



VI. EXCHANGE IN ORDERED SPIN 
BACKGROUND 

It is to be noticed that unlike our previous SCR for- 
mulation^'i the action no longer has full spin-rotational 
invariance, owing to the anisotropy. (Notice that putting 
J2 = recovers the familiar rotationally invariant 
Anderson-Hasegawa form.) Hence, even if we start with 
a perfectly ordered spin background with all spins par- 
allel, we must also choose an orientation relative to the 
anisotropy axis. We consider 3 cases. In the following, we 
replace epe = 0, and tMo = A for convenience. Hence- 
forth, J2 will mean z J2, where z is the coordination num- 
ber. 

First, we consider the case 6i — {) for all sites. Then, 
translational invariance is restored, and we can write the 
relevant terms of the action [S] in momentum space and 



+Ek-nF{Ek-) - (A + J2)iif{A + J2)] e 



k-{fi—rj ) 



(9) 



where Ek± = ^+-h±y/iA+j2) + ^^^ This gives the final 
expression for the exchange which has to be evaluated on 
a square lattice. 

The case of 9i = tt gives, proceeding in a similar fash- 
ion, 



k+) 



+E_nF{Ek-) - (A -f J2)nF(A + J2)] e 



(10) 



The same as = 0! This is not surprising, since uniaxial 
anisotropy does not distinguish between 6 — and 9 — n. 
Now, let us consider the case of 9i — ^. 



kn 

Hence, 

1 



IUJ„ 



1 

^ '-FM 



1 



2 iujn — A 2 iuJri 



u ^ 



^ 2 

kn 



el{2iuj„-A) el{2iuj„-A) 



{iLOn - Af 



Af 



iLlkn (11) 



7L7fcn)l2) 



Thus J2 drops out entirely! Thus, the expression for the 
exchange here would coincide with that for J2 = eval- 
uated earlier^-, which is also not surprising because the 
anisotropy axis has no component along the perpendicu- 
lar direction, so that the superexchange has no effect on 
spin configurations along this direction. Replacing the 
Green's function. 



U 



'(2iu;„- A) 
_{iujn-Af 



(13) 
(14) 
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This upon evaluation of the Matsubara 
sum gives the famihar expressioni^ U — 
J2k [Eu+nF{Ek+) + Ek-npiEk-) - AnpiA)] as be- 
fore. This is to be expected, since the anisotropy, and 
consequently, the superexchange, will have no effect in a 



direction perpendicular to the easy axis. 



Finally we write the exchange for a general 6i = 9. 



''FM 



cos'^^(2iuj„ 



J2COS6) sirr^{2iujn — A + J2COS9) 



{iujn — A — J2COS6 



{iuJn — A + J2COS6 



X Gkn 



(15) 



where 



el[{iuJn- A) + J2COS^0] 

(iw„ - A)2 - Jlcos^e 



(16) 



VII. RESULTS 




The results of the exchange calculations are shown 
here. For the case of 9 = (up) the exchange is as 
shown in Fig[TJ Since we are only interested in looking at 
the ferromagnetic phase, hence only the dominant near- 
est neighbour exchanges are shown. As explained earlier, 
the same graph is obtained for the 9 = tt (down) case. 
For comparison, the J — curve for the same A is also 
shown^'^. This may also be thought to be the exhcnage 
data for the 9 ^ tt/2 case. Two effects are immediately 
observed. Firstly, the extent of the ferromagnetic phase 
in the filling regime increases for the finite superexchange 
case, as compared to no superexchange. This is proba- 
bly why for an extended filling regime, the ferromagnetic 
behaviour persists rather than the antiferromagnetic be- 
haviour, in the 3d-5d compounds. In particular, while the 
exchange already becomes positive for a filling of n w 0.7 
for the J2 = case, it continues to be positive way be- 
yond n — 1 for finite J2 ■ Secondly, the magnitude of the 
negative part is also larger, proving that the ferromag- 
netic Tc is enhanced by superexchange. Another inter- 
esting point to note is that the exchange continues to be 
large and negative till about n = 1, which corresponds to 
iV = 3 in the real double perovskites owing to t2g degen- 
eracy (eg. Sr2CrOs06), while its magnitude diminishes 
beyond that. The point n = 2 shows a strong isolated 
antiferromagnetic tendency, possibly due to the filling of 
the Mo level with opposite spin in the half-metallic state. 

Next, to bring these results in perspective, and estab- 
lish them on a sound footing, we also performed extensive 
numerical simulations of the original spin-Fermion model 
Hamiltonian, using the technique of exact diagonaliza- 
tion coupled with Monte Carlo (ED-I-MC). Calculations 
performed in real space on a 8 x 8 cluster is shown in 
the figure Fig [51 We used A = —2 and J2 = 1 as before. 
It is found that inclusion of the J2 term in the Hamil- 



0.5 1 1.5 2 2,5 3 

Figure 1: Exchange vs filling obtained from SCR, compared 
between cases with and without superexchange 



tonian not only results in a tremendous increase in the 
Curie temperatures, but also in the proliferation of the 
ferromagnetic phase to higher fillings. In fact, for these 
regime of parameters, the ferromagnetic phase is found to 
make substantial entries into those regions of filling which 
were earlier reserved for the antiferromagnetic phases^. 
However, both effects are far more augmented than sug- 
gested by the SCR calculations. Of course, the SCR is 
performed assuming ordered spin backgrounds appropri- 
ate at low temperatures and also only one loop of the 
actual SCR process is performed to derive the analyti- 
cal expressions given in this paper, hence it should only 
be considered as a pointer to the actual renormalization 
effects close to Tc- Interestingly, in the exact numerical 
simulations, both the ferromagnetic lobes at low and high 
filling join together in the case of finite J2, giving a max- 
imum Tc at N=l (which corresponds to N=3 considering 
the t2g degeneracy in the real material). This is con- 
sistent with the enigmatic high Tc observed in the com- 
pound Sr2CrOs06-^-'^^. Also, it substantiates the view 
that this compound is at the threshold of a magnetic 
and electronic transitiort^^. Thus the mechanism of su- 
perexchange coupling between B and sites provides 
an alternate, and more acceptable mechanism for the Tc 
increase, rather than the Hubbard U mechanism invoked 
earlier—. 
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A=-2 




Figure 2: Tc vs N ferromagnetic part of pliase diagram from 
8x8 cluster compared between J2=0 and J2 = l (A = —2) 



microscopic Hamiltonian, we showed that the superex- 
change is the main player in enhancing Tc as weU as 
stabihzing ferromagnetism in the Cr-compounds. Spin- 
orbit couphng breaks the rotational symmetry on the B 
site, so that different orientations of spin configurations 
have to be considered for the SCR. Amongst the configu- 
rations considered, those along the easy axis are found to 
increase the effective nearest neighbour ferromagnetic ex- 
change, as well as contribute to the increase in the filling 
extent of the ferromagnetic phase. Both are signatures 
of stabilization of the ferromagnetic phase, which results 
in the increase in Tc, and continuation of the ferro phase 
to higher fillings, as for Rh and Os compounds. Exten- 
sive numerical simulations are found to reproduce both of 
this behaviour, and provides a sound backing to the pro- 
posed model Hamiltonian as an useful one for describing 
the behaviour of Cr-based 3d-5d double perovskites. 
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Wc have deliberately considered a larger value of J-z than 
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tual values of J2 may leave an window for the antiferro- 
magnetic phases as well. The effects of phase competition 
in this Hamiltonian will be considered in a later work. 



